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Abstract 

This paper provides a discrete time LIBOR analogue, which can be used for arbitrage-free 
discretization of Levy LIBOR models or discrete approximation of continuous time LIBOR 
market models. Using the work of Eberlein and Ozkan _3j as an inspiration, we build a 
discrete forward LIBOR market model by starting with a discrete exponential martingale. 
We take this pure jump process and calculate the appropriate measure change between the 
forward measures. 

Next we prove weak convergence of the discrete analogue to the continuous time LIBOR 
model, provided the driving process converges weakly to the continuous time one and the 
driving processes are PII's. 

We also demonstrate the results of a practical implementation of this model and compare 
them to an implementation of an arbitrage free discretization by Glasserman and Zhao [J]. 

1 Introduction 

LIBOR Market Models are intensely used in banking and finance since their appearance in the 
paper by Brace, Gatarek and Musiela [1 . A rich literature has developed in the field, especially 
the book by Brigo and Mercuric 2 should be mentioned as a accessible reference. Since all prac- 
tical calculations need discretization of some sort, arbitrage free discretization as in Glasserman 
and Zhao 4 is desirable. With the papers by Jamshidian [6 and Eberlein and Ozkan [3] more 
general driving processes became the center of attention. The dissertation of Kluge J, offers 
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an insight into application and calibration of those models. However in continuous time, one is 
usually forced to carry out the so called " frozen drift approximation" . In a time discrete model 
driven by discrete random variables, such an approximation could possibly be avoided, which 
was the original motivation to construct such a model. 

It thus seemed natural to look for a fairly general discrete time analogue of these models. It 
should be arbitrage free and the analogy to the continuous time case should be well apparent. 
We treat this in section Q. 

Then we devote an equally important section to the convergence of a sequence of discrete ana- 
logues, under certain integrability and convergence assumptions, to a continuous time model. 
This property has been the starting point of our interest in discrete LIBOR models. We get a 
convergence result and can demonstrate, that our discrete model can be applied to find approx- 
imations for the exact (non frozen drift) Levy LIBOR model. 

The next section ([2]) gives an overview on the results from a model implementation and compares 
the results to a model based on the work [4 . This is carried out in more detail in section ([5]). 
The final section is devoted to an outlook on further possibilities for work based on these results. 
All proofs can be found in the appendix 

Apart from the papers above, the methods we use can be found mostly in the books by Protter 
|10) and Jacod and Shirjaev The theorems relevant for the convergence result can be found 
in appendix (B\ . 

If the readers interest is purely in the implementation results, I recommend to look at sections 
3.2, 4.1 and 5. If the readers interest is more in the convergence result, sections 3 and 4 cover 
this entirely. 

2 The Merits Of The Approach 

We can show that there is a discrete LIBOR market model analogue such that all LIBOR rates 
converge under their respective forward measures and also jointly under the terminal measure 
to a continuous time model, as soon as the driving process converges. 

The implementation of such a model gains us the possibility to calculate without the frozen drift 
approximation, assuming the driving process is simple enough (see section (O for details). 
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Below you can see the results of 2 implementations of this approach, compared to an implemen- 
tation of one arbitrage free discretization by Glasserman and Zhao from and implied volatility 
smiles derived from the data. The 3 models have been chosen, such that all 3 will be arbitrage 
free and converge to the lognormal LIBOR market model, without frozen drift. 



LogBernoulli Discrete Model Implied Volatility Lognormal Discrete Model Implied Volatility Glasserman -Zhao Model Implied Volatility 




If we look at the results from left to right, model 1 (exact calculation with non frozen drift) 
yields considerably higher volatilities than the other 2 models. Model 2 (log normal discrete) 
yields higher implied volatilities than the Glasserman-Zhao discretization for caplets with small 
strikes but smaller volatilities as the strike grows. Model 3 is one of the the Glasserman-Zhao 
arbitrage free discretizations from [1]. 

The details on the implemented models and the convergence result for a refined time grid can 
be found in the following sections. 

3 Model Approach And Properties 
3.1 Aims 

Based on what goals did we derive the discrete time LIBOR market analogue? 

• We expect the model to fit a given finite tenor structure with initial data for each point in 
the tenor structure. 

• We expect the model to give an explicit expression for the evolution of the respective rates 
in the future. 

• We expect the model to prevent arbitrage opportunities between different LIBOR rates. 
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• Finally, wc want the model to approximate continuous time LIBOR market models if we 
let the discrete time grid become infinitely fine. 

All of this can be accomplished using the following model: 



3.2 Discrete LIBOR Market Model 

We assume as given a stochastic basis {ft, F, P) with F = {J-'t)t>o satisfying the usual condi- 
tions. 

Given a time horizon T*, we represent running time as a finite grid of m = {n+l)p positive real 
numbers {n,p G N) 

ti:=—T*, i = 0,l,...,m 
m 

and we define an (n + l)-element tenor structure through 

T,:=^T*. i = l,...,n + l. 

We further assume we have a family of positive real numbers {Xij)i=o,...,m,j=i,...,n and positive 
real values {Z/(0,T,)}"^o- Furthermore let 5 := 

Definition 1 (Discrete Forward LIBOR Market Model) A family of stochastic processes 
{I/(fi,Tj)}i=o,...,pj,j=i,...,n o,nd a family of equivalent probability measures {Pj}j=2,...,n+i consti- 
tutes a LIBOR market model if 

1. L{., Tj) is a Pj+i martingale for j = 1, . . . ,n. 

2. For j = 2,...,n 

dFj ^ _ l + SL{t^,Tj) 



d¥j+i^^'^ ~ l + SL{0,Tj)- ^^-^^ 
Indeed such a model can be constructed by 

Theorem 1 (Existence Of A Discrete LIBOR Market Model) Given an adapted time dis- 
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Crete process build from random variables (^tj)i=i,...,m such that 

m n 

i=i j=i 

There exists a discrete LIBOR market model (i(.,Tj))"^j 

i 

L{ti, Tj) = L{0, Tj) exp A„,(Xt„ + 6^J), i = l,...,jp (3.3) 



= -T^ logEp,^! (exp(Ai,XtJ|jr,^_ J. (3.4) 

The proof is straight forward. It can be found explicitly in the appendix. 
This approach allows us to define implied Bond price and Forward dynamics: 
Let viU) := u for ^T* < ti < ^T*. 

Definition 2 (Bondprices And Forward Price Process) We define the family of processes 
called bond price processes (-B(ti, rj))i=o,...,pj,j=i,...,n from the LIBOR rate processes ^I/(.,Tj)^ 
through 

and then 



iJ(i„T,)=iJ(0,T^(,,)_i) n i + 5L(tM - Vti^T,,r = l,...,n (3.6) 

u—r]{ti) — l 

We call 

FB{ti,Tj,Tj+,):=l + SL{ti,Tj), Tj G {Ti, . . . ,T*},5 = (3.7) 
the forward price process. FB{ti,T,T + 6) is aVx+s martingale. 
With that, we get a full LIBOR Market Model. 
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4 Convergence Result To Continuous Time LIBOR 

While we view the discrete LIBOR market model analogue to be of intrinsic interest, the central 
idea when developing the model was to allow for an approximation of continuous time LIBOR 
market models without having to freeze the drift. 

In this section we therefore derive a convergence result of a sequence of discrete time models 
to a continuous time model. First, we derive a difference equation representation to make the 
applicability of convergence results more immediate in the following. 

4.1 Difference Equation Presentation of the Dynamics 

We apply the following definitions to the process Xj. + Bt^ = X]l=i(^tu + ^t^)- 
Definition 3 (Jump Measure) We call 

fiiA, {tjti) := e {U}to\Xu + K € A} j>0 (4.1) 
the jump measure of the process Xf . + Bj . . 

With this definition it holds for / continuous and E[f{Xt-)] < oo,Vi 



r [ /(x)^(dx, dt) = V f{Xt^ + bt^). (4.2) 
Jo Js. 



Furthermore we need 

Definition 4 (Compensator) We define the compensator of a random jump measure ji with 
respect to a measure P to he 

KA,{i.K=i) =Ep(MA,{i.K=i)l-^*.-J (4.3) 

One can then show 

Theorem 2 (Dynamics And Measure Change) Given a discrete LIBOR market model as 
constructed in theorem ([T]). Then the dynamics of the LIBOR rates {L{ti,Tj)\ 

V / 2— 0,...,jp,j — l,...,n 
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under their respective forwardmeasures (Pj)j=2,...,n+i o-fs- 

AL{ti,Tj) = L{ti-i,Tj) / (exp(Aija;) - - u^+^){dx,ti) under P^+i (4.4) 
Js. 



where 



and 



{ti}U) ■■= *{U e {ti}U \Xu + K e A} (4.5) 



u=+\A, {ti}U) = Ep,+,(m(A {ii}r=i)l-^t.-J- (4-6) 
And then the measure change from Pj+i to Pj is reflected in a change of the compensator 

u^idx, {ti}) = ie{ti-i,Tj){exp{Xij{x)) - 1) + iy+^{dx, {ti}). (4.7) 

Here i is completely analogous to continuous time: 

Due to the difference equation representation, we can write the measure-change down as a change 
in the compensator, as in the continuous time case. 

4.2 Approximation Theorem For Discrete LIBOR 

For the convergence result of the sequence of discrete time models, we embed the discrete models 
in a common space. This is best describes the way we think about the discrete time models if 
they are used for approximation of continuous time models. 

Theorem 3 (Approximation Theorem) Given a sequence {X^'^^)keN of discrete time PIIs 
on stochastic bases B^''^ = (fiC^), ,7"^, (J'f ^)»6N,Pi+i) • 
For each k we define an adopted sequence of random variables 

C/f := ^ - X.^^i, i G N (4.9) 

and a change of time (a*) to consider then as processes with paths in D([0, T*]) := {/ : [0, T*] — )• 
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i\fis cadlag} with fixed jump times the processes 

X«:= ^ ^/^ te 



+ ■ 

l<i<o-f 



As k goes to infinity we let the mesh of the set of jump times go to 0. 
Now we define a PII (^X^ on (J't)t, P„+ij 

Xt + Bt^ I bsds+ I I cldWs+ I I x{fi-iy){dx,ds), (4.10) 
Jo Jo Jo Jo Jm 

where {J^t)telo,T*] is the filtration generated by (^t)te[o,T'] OL'^'d 

Bt^ f b,ds. (4.11) 







Also (Xf) has to fulfill the conditions on a driving process for a Levy LIBOR model as in 
This imposes the following condition on the drift characteristic of the driving Levy process 



Xis,T* - 6)bsds ^ f Cslx{s,T* -Sfds+ f f (e^'^''-'^'^^^'' - 1 - X{s,T* - 5)x)v{ds,dx). 

Jo 2 Jo "'R 

(4.12) 

We assume that 

XC") + ^(fe) X + B, weakly as processes (4-13) 



for 



B^=Y.^ (4.14) 

i=i 



(where the b'^^^ are defined m l3.4|) . and it holds that 

sup E(/ / exp{ux)F^''\dx)ds) < oo,u< {1 + €)M,M >y2\K-^Tj)\, (4.15) 

fcSNUoo Jo J|a;|>l j^-^ 

v^^\dx,ds) = Fl^\dx)ds.,v°°{dx,ds) = F^{dx)ds = v{dx,ds) = Fs{dx)ds. 
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as well as 

sup E{[ [ {x^ A 1) F^'''>{dx)ds) <oo. (4.16) 

We define discrete LIBOR market models {L{ti, Tj)'--''^)i^Q i jp'' j=i,...,n.keti> driven by the (X''''))^^^ 
as in the section on discrete forward LIBOR modeling. In other words the rate L{., Tj)^''^ is driven 
by the exponential of (X*^'^^) minus a drift, making the exponential of (X^'^^) into a martingale 
under p'^i . 

We assume that for := sup^^g i jpk < ti it holds that A^j.^. X{ti,Tj) pointwise, and 
X{.,Tj) : — ^ are linearly bounded, locally Lipschitz and positive functions. Furthermore 
we assume we are given starting values lL{0,Tj)] which are strictly increasing in j then 

V / j—l,...,n 

L{.,Tj)^''^ under Pj^i H-^Tj) under Pj+i VI < j < n (4.17) 
weakly as a process. 

Lemma 1 (Contiguity of forward measures) Under the conditions of the preceding theorem 
any two sequences of forward measures are contiguous. 

5 Implementation Results 

In order to support the practical use of the discrete time LIBOR market model, we calculated 
caplet prices and volatility smiles for a reference rate. This shows that not only is pricing of 
options feasible based on the discrete time model, but the results allow to find other arbitrage 
free prices for discrete approximation of the lognormal LIBOR market model. 
The discrete model ([T]) was implemented using 

1. A purely discrete driving process with Bernoulli distributed random variables (Xt-) with 
parameter p. For this model, it was well possible to calculate the exact LIBOR model. 

2. A model driven by standard normal (Xt^) . This model was used for Monte Carlo simula- 
tion, as the exact local drift would be to costly to calculate. 

As a reference, an implementation of one of the models described in using the same volatili- 
ties (Ay ), was carried out. 
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5.1 General Setup 

Our aim was to calculate the price of forward caplets on LiTc^^T^) in a model with 10 LIBOR 
rates. Hence the tenor structure was 

T,^i. i = l,...,ll (5.1) 

with the Ti denoting years and the step size being 6 = 1 giving us a time grid t = 0, 1, 2, . . . , 11. 
The time Tn was not used to model a LIBOR rate but only to define the terminal measure 
and hence the forward measure for rio))tiGO,....io- 
Bond prices B{Q, Ta), . . . , B(0, Tn) were calculated from ((3?6)) . 
The caplet price to be calculated was 

n = 5(0, r6)Ep, [i(T5, Ts) - KL{Q, Tg)] (5.2) 



for K values 





K2 






^5 








0.6 


1 


1.4 


1.8 


2.2 


2.6 


3 


3.4 



The initial data was 



i(o,ri) 


i(0,T2) 


i(o,r3) 


L(0,T4) 


^o,r5) 


i(0,T6) 


L(o,r7) 


i(o,r8) 


L(0,T9) 


i(o,rio) 


0.0207 


0.23 


0.0262 


0.28 


0.0292 


0.0318 


0.0342 


0.0362 


0.0379 


0.04 



For comparison purposes, all chosen models have driving process variables X(. with E[X4j = 

and V[XtJ = 1. All models converge to the lognormal LIBOR market model for refined time 
grids. 



As a measure of how the model results related to one another, we compared the caplet prices, in 
terms of implied volatilities, obtained for the below a priori volatility structure for all 3 models. 
The a priori volatility parameters were 



Volatilities {Xij) 


Ail 


Ai2 


Ai3 


Ai4 


Ai5 


Ai6 


A.7 


Ais 


Ai9 


Aiio 




0.34 


0.32 


0.3 


0.28 


0.26 


0.24 


0.22 


0.2 


0.18 


0.16 
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The volatility parameters were independent of the time index i and depended on the rate index j. 

5.2 Implementation Of The Log-Bernoulli Model 

Based on the central equation from section ([3]), Theorem ([Ij, equation (j3.3|) 

■i 

Liu, T,) = L(0, Tj) exp K,{Xt^ + 6^J) (5.3) 

u=l 

using Bernoulli Variables (^tj, IP[-Vt. = 1] = 0.5 and ^[Xt^ = — 1] = 0.5, it was possible to 

calculate all paths of the resulting discrete model. The pricing of a caplet then reduces to 

summing up over all possible payoff values times the respective path probability. 

While the number of paths is 2^, in this example, and we have to calculate from the terminal 

rate to the target rate L(, 5), the speed of the computation is fast for the sample data. 

The main computational issue is to properly calculate 

K = "T^ logEp^^,(exp(A,jXtJ|J-t^_J = (5.4) 

-^logIEp„,,(cxp(A,X,) n 
Differently put, we need the compensator under the terminal measure P„+i 

n 

v^{dx,{U}) = (J](^(t,_i,rfc)(exp(A,fc(a;)) - 1) + l)W+\dx , {U}) . (5.5) 

We therefore need a different compensator value depending on time and the £{ti-i,Tk) of the 
earlier rates. For a model with 10 rates this proved to be no issue, once the compensator product 
was implemented. Our implementation is based on the fact that 

n 

([](£(t,_i,Tfc)(exp(A,fc(x)) -l) + l)))i."+i(dx,{tj) = (5.6) 

k=j 
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where a denotes all subsets of the set {j, j + 1, . . . , n} and ct*^ is the complement of a in {j, j 
1, . . . 

The resulting implied volatilities are 



K 


0.6 


1 


1.4 


1.8 


2.2 


2.6 


3 


3.4 


implied Volatility 


0.542 


0.396 


0.334 


0.32 


0.277 


0.2875 


0.2832 


0.25 



In the figure below, the huplied volatihties are rescaled by a factor of 100. 

Volatihty Smile for the exact discrete time log-Bernoulli model 

LogBernoulli Discrete Model Implied Volatility 



1.0 1.5 2.0 2.5 3.0 3.5 

Strike as factor (limes initial Value) 



5.3 Implementation Of The Lognormal Discrete Analogue 

We again used section ([3]) , theorem ([TJ , equation 



L{U, T,) - L{0, T,) exp K,{Xt^ + 6^J) 



M = l 



(5.7) 



but here with Xt^ being standard normally distributed. The compensator under the terminal 
measure is once more 



v^{dx,{U}) = C[{{l{U^^,Tu){eM^^k{x)) - 1) + l)))y^+\dx,{U}) 

k=j 



(5.8) 
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but here £{ti^i,Tk) can have infinitely many possible values and the path by path calculation 
of the above model becomes impossible. Therefore, pricing was carried out, using 500000 paths 
from Monte Carlo Simulations. This model is somewhat slower than the Glasserman-Zhao dis- 
cretization in the present implementation, due to the compensator product. It however yields 
different results and therefore yields other arbitrage free prices. The resulting implied volatilities 
are 



K 


0.6 


1 


1.4 


1.8 


2.2 


2.6 


3 


3.4 


implied Volatility 


0.518 


0.333 


0.276 


0.253 


0.241 


0.234 


0.23 


0.226 



In the figure below, the huplied volatihties are rescaled by a factor of 100. 



Volatility Smile for the lognormal disrete model - Monte Carlo Simulations 

Lognormal Discrete Model Implied Volatility 



Strike as factor (times initial Value) 



5.4 Implementation Of One Glasserman-Zhao Arbitrage Free Discretiza- 
tion 

For the arbitrage free discretization of [3], p. 43 we used the equations (21), (24) and (25) with 
6 = 1 and h = 1 
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with the same i,j as in the above models except for (L(ti, riO))t.go,....io which is modeled arbi- 
trage free already. Here 



= W,{i - 1) exp ( - \^a]{i) + a,{i)Y, 



(5.10) 



with 



crj(i) = Ay + XI Y" 



(5.11) 



and A/'(0, 1) and the (li)i=i^...^io are independent. 

The pricing formula we used was from on page 45 adapted for iV = 10 and n = 5. 
Here Oj(i) again includes local terms, dependent on the value for a higher maturity W^+i. 
Therefore, for a reference implementation, we also used Monte Carlo Simulation for pricing 
(500000 paths). There is no compensator product here and compared to the log normal discrete 
model, this method is faster. The log-Bernoulli model however turns out to be even faster (even 
using exact calculations) than this method. Test for further data will follow in future work. 



K 


0.6 


1 


1.4 


1.8 


2.2 


2.6 


3 


3.4 


implied Volatility 


0.437 


0.313 


0.285 


0.2745 


0.269 


0.265 


0.263 


0.261 



In the figure below, the implied volatilities are rescaled by a factor of 100. 



Volatility Smile for one arbitrage free discretization of Glasserman and Zhao 

Glasserman-Zhao Model Implied Volatility 




Strike as factor (times initial Value) 
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6 Summary And Outlook 

We have proven the existence of the arbitrage free discrete LIBOR analogue ([Ij and its con- 
vergence to a continuous time model under reasonable assumptions (convergence of the driving 
processes, uniform integrability) . 

Further work will be along the lines of the work by Musiela and Rutkowski (S). We aim to 
establish a complete discrete LIBOR framework, namely forward models and the implied savings 
account. Also we want to adapt the convergence result (jSj for joint convergence under terminal 
measure. 

Establishing some relation to the general HJM framework would be exciting. 

Interest in LIBOR Market Models is high as ever as the work of Keller-Ressel, Papapantoleon 

and Teichmann ^ and recent work by Eberlein et al show. 
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A Proofs 

Proof of theorem [H section ([3]): 

First we check what L(., Tj) has to fulfill to be a Pj+i martingale. 

We proceed inductively: Assume ( L(ti, Tj) ] has been constructed to be a martingale. 

V / i=0,l,...,u-l 

We look at L{tu,Tj) 

IEp,-+i(i(in,r,)|J-t„_J - i(t„-i,Tj)Ep^.+i(exp(A„,(Xt„ +6^J)|J-t„_J. (A.l) 

We therefore have as martingale condition 

Ep,+i(exp(A,y + blJ)\Tt,^_J = 1 o Ep^._^,(exp(A„jXf JlJ-i^i) = exp(-A„j&{j. (A.2) 

Here we choose to be ^-measureable. 
This yields 

bl = -^logEp^+,(exp(A,jXiJ|J-t_J. (A.3) 

So we have that L{tu,Tj) fulfills the martingale condition with this choice of drift. 

Now assume we are given {L{ti,Tj})i=o,...,nj under Pj+i as a Pj+i martingale. We define Pj 

through 

dPj _ l + 5L{U,Tj) 



dPj+i' ' l + SL{0,Tj) 

All appearing expectations exist due to the integrability condition (13. 2p . □ 
Proof of theorem section (H]): 



We have that 



AL(t„T,-) =i(0,Tj)exp(^A„j(Xt„ +6^j)(exp(A,,(Xt^ +6^J) -1) = (A.5) 

M=l 

L(i,_i,r_,-)(exp(A,,(Xt^ + blj) - 1). 
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And with the definition of /i and i^^^^ it holds that 



(exp (A,, (a;)) - l)fi{dx, {U}) = exp (A,, (X^, + 6^) - 1 (A.6) 



and 

(exp(A,j(a;)) - iy+^{dx,{U}) = -exp(O) + 1 = 0. (A.7) 



So 

ALiU,T,) = L{U^uT,) / (exp(A„-(a;)) - - v^+^){dx,{t,}). (A., 



From the definition of v^^^ it follows that 

{e^^{\,x)u'+\dx,{U}) = 1. (A.9) 



We then proceed to look at the forward price process F{ti,Tj,Tjj^.i) = 1 + dL{ti,Tj). The 
dynamics under Pj+i are derived as 

AFiU,T,,T,+i) = 5ALiU,Tj). (A.IO) 

We have described AL{ti,Tj) above. Therefore we get an analogy to the SDE of the forward 
price in continuous time 

AF{U,Tj,T,+i) = F{U^i,Tj,Tj+i) [ e{U^i,Tj){exp{X,j{x) - - iy^+^){dx,{U}). (A.U) 

JR 

From the proof of theorem ([1]) we know that 



dFj+i'" Fb{O.Tj,T,+i) 

With the discrete Girsanov Theorem and the stochastic difference equation representation of 
Fb [U , Tj , Tj+i ) , we get 

'^^ -{t^ = Zi=£{Mi) (A.13) 
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with 



Ml^y^f ^(t„_i,r,)(e^-(-) - l){^,-u^+'){dx,{U}). (A.14) 
We apply discrete Girsanov theorems to get 

Y^{Xi) = / (^(t,_i,T,)(e^-(-) _ 1) + _ uO+^){dx,{U}). (A.15) 



Therefore 

yi{xu x^-i;xi) ^ l{U^i,Tj){e^^=^''^^ + (A.16) 

We have that 

dJlixi,. . . ,Xi-i]Xi) ^ v'{dx,{ti}) A dff^^{xi,...,Xi-i\Xi)^v'^'^{dx,{ti}) (A.17) 

and 

df^ 

yi{xi, . . .,Xi-i;x^) = j+i {xi, ■ ■ .,Xi-i;xi). (A. 18) 

This gives us 

i^^dx,{U}) = (^(i,_i,r,)(exp(Ay(x)) - 1) + iy+\dx,{U}). (A.19) 

And this concludes the proof. □ 
Proof of theorem section Q: 

We denote by 

:= inf e R+ [af = i} (A.20) 

the embedded grid points of the discrete time models. Taking limits from the righthand side, we 
get weak convergence everywhere, after we have proved the result for the grid points. 
With condition 

sup E( / / exp(ua;)Fi'=)(da;)ds) <oo,u< (l + e)Af,M> V|A(.,T,)|, (A.21) 

fceNU+oo Jo J\x\>l 



1=1 
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fulfilled, there exists a continuous time LIBOR model L{t,Tn) with driving process X as in ^ 
and discrete Models {L{ti,Tn)^'''^)keN with increments U^'^''\ 

The functions A(.,Tj) are linearly bounded and locally Lipschitz and the (A|^-')j^Q p(k}j „ 
are pointwise convergent, deterministic functions, which we need for ^ in Jacod and Shirjaev 
([5], page 578, theorem 6.9). We have the stochastic difference equation. 



AL«(t„T„) =iW(t,_i,T„) / (exp(A,„(a; + 6r/'^)) - - 'y"+'^'^''^){dx,U). (A.22) 



To see that this converges to the SDE of L{t,Tn) we need conditions on the convergence of 
stochastic differential equations in a general semimartingale setting. Such conditions can be 
found for instance in Jacod, Shiryaev [S], page 577f . 
First we need certain existence and uniqueness results: 

In our concrete case and using the notation from the appendix, we have that: 




Y = Z + f{Y^)-0 



with 



^(fe) ^ _^(fc)(0,T„) VyfceN 




Vfc e N 




and 



• Our O^^^ is weakly convergent by assumption to 



O := 



/(exp(A(i,r„+i)a:) - - u'-+^){dx,dt) 



(A.23) 



Therefore certainly tight. In absence of a drift component, we automatically have that 
the drifts variation norm V^ar(B"'(''))t (see [5], page 27, 3a.) is tight. Therefore {0^^'^)k is 
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P — UT (see 0, page 380, 6.15 and page 382, 6.21 or appendix B, definition ([SJ). For the 
definition of P — UT see [S], page 377, 6.1 or appendix B, definition ([5]). 

• We have that /'^''^ is Lipschitz and linearly bounded for all (k). 

• We have assumed O^'^^ — O weakly and shown that it is also P — UT and know that 
ZC^) = L{0,Tn) L{0,Tn). Therefore we have (O^*^) , Z^^)) ^ {0,Z) weakly and by the 
theorem on weak convergence of solutions of SDE's we have (0('=),zW,rW) ^ {0,Z,Y) 
weakly. 

Those conditions fulfilled entail that the sequence of discrete processes converges weakly to a 
weak solution of the proper SDE in continuous time. 
That SDE is 

dL(i,r„) = L{t^,Tr,){cldWs + [ (e^(*'^")- - l)Oi - (A.24) 



By the same reasoning, we get that AFB{ti,Tn,Tn+i) also converges to the SDE of the forward 
in continuous time. 

By contiguity of the forward measures (see section ([4]), lemma ([ij and the proof below) all further 
rates converge under their respective forward measures. □ 
Proof of lemma ((!]), section (|3]): 

We denote by 

t'^ := inf e R+ = i} (A.25) 

the embedded grid points of the discrete models and prove the result for the grid points. 
The measure change between two consecutive forward measures is given by 



'^i (f\ _ _ F(ATi\ - ^B{t,Tj,Tj+i) 

^(t) -Z,- £[M, ) - ^^^^^^^^ 



and with respect to the terminal measure 

<^3 _fr ^s(i,Ti,r,+i) 



(A.26) 



^ — 3 
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Let Z^'^''^ denote the measure change from P^*j_\ to Pi+i- 

With this notation we use [5] p. 288, theorem 1.10, (m'). Our sequence is uniformly integrable 
due to (I4.15|) . This yields contiguity according to [5 p. 288, lemma 1.10, (m'). □ 



B Weak convergence of processes 

Definition 5 (Predictably Uniformly Tight) For each integer n let (f^", J"", F", P") 

be a stochastic basis. We denote by % the set of all predictable processes iJ" on i3" having the 
form 

k 

i=l 

with fc G N, = So < si < ... < Sfc < Sfc+i and Yi is J- ".-measurable with \Yi\ < 1. 
We can define an elementary stochastic integral by 

k 
i=l 

Now a sequence (X") of adaptedf on their respective Bases) cadlag d-dimensional processes is 
P-UT if for every t > the family of random variables {J2i<i<d-^^'^ ' • ^ N,i/"'* G H") 
is tight in R meaning 

d 




> a) = 0. 



Assume the following setting: 



• We are given an equation 

Y^Z + f{Y_)-0 

as above 

• For each n we have a stochastic basis = (il", J^", F" , P") and an equation 



y" = z" + /„(y")-o' 



(B.l) 
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• O" is a d-dimensional-semimartingale on that basis. 

• Z" is a q-dimensional cadlag adapted process. 

• /„ are functions M'? ^ M'' x M'^ such that each equation above admits a unique solution. 
Then there holds 

Theorem 4 (Weak Approximation) Assume the functions /„ fulfill Lipschitz and linear bound- 
edness with constants not dependent on n and /„ f at least pointwise. Assume further that 
the sequence O" is P-UT. Then ifY"' denotes the unique solution of the sequence of equations 
there holds 

If{0'^,Z'') {0,Z) weakly, then (0",Z",y") ^ iO,Z,Y) weakly. 
Proof: 

See Jacod and Shirjaev [5], page 578, theorem 6.9. □ 

In the case of the theorem on discrete approximation we need to show especially P — UT of 
the driving process of the SDE. However we have weak convergence of the sequence of driving 
processes already and in that case there holds 

Theorem 5 (Tightness And P-UT) Let (0")„ be a sequence of d-dimensional semimartin- 

gales with characteristics and second modified characteristics (_B",C",t^") and C". 

// the sequence (0")„ is tight, then the sequence is P — UT iff the sequence Var{B^'^)t is tight. 

Proof: 

Using Jacod and Shirjaevs book [5] we combine page 380, 6.15 and page 382, 6.21 and the theo- 
rem follows. □ 
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